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Preface 

The  purpose  of  this  thesis  was  to  analytically  predict 
and  experimentally  measure  the  flow-induced  pressure  field 
in  a  reed-valve  arrangement.  Reed  valves  are  typically  used 
in  sealed  compressor  units.  The  thrust  of  this  work  was  to 
find  an  analytic  solution,  compatible  with  a  small  computer 
system,  for  a  "first-look"  at  the  pressure  field  on  the  up- 
flow  side  of  the  reed  valve.  Analytically,  the  flow  was 
modeled  as  one-dimensional,  piece-vise  isentropic  with  the 
possibility  of  an  imbedded  normal  shock.  The  experimental 
portion  of  the  thesis  work  involved  operating  an  idealized 
reed-valve  model  at  various  steady-state  pressures.  The 
pressure  field  beneath  the  valve  plate  was  measured  using  a 
computerized  data  acquisition  system  which  relied  on  a 
scanivalve  as  the  primary  pressure  gauge. 

An  interesting  side  note  identified  during  this 
research  was  that  if  you  make  a  simple  model  of  the 
conf iguration,  using  a  spool  of  thread,  a  pin,  and  a  piece 
of  paper,  you  can  create  a  supersonic  flow  and  a  normal 
snock  by  blowing  through  the  spool. 

Several  people  gave  very  valuable  assistance  thoughout 
this  project.  My  advisor,  Lt  Col  Eric  Jumper,  provided 
insight  into  the  analytic  solution,  provided  direction,  and 
shared  in  the  enthusiasm  of  the  experimental  work.  Nick 
Yardich  and  Leroy  Cannon  were  ever  ready  with  their 


professional  advise  and  assistance. 


A  special  thanks  is  due 


to  Jay  Anderson  who  diligently  assisted  me  in  setting  up  the 
computerized  data  acquisition  system.  The  members  o±  my 
class,  GAE-B5D,  were  very  supportive  and  I'm  proud  to  call 
each  one  of  them  a  friend.  finally,  my  love  and  thanks  are 
devoutly  extended  to  my  wife,  Alice.  Her  encouragement, 
professional  critique  and  support  were  invaluable  even  when 
the  Atlantic  Ocean  spanned  the  distance  between  us. 
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Abstract 

The  flow-induced  pressure  field  in  a  reed-valve 
arrangement  was  analytically  predicted  and  experimentally 
measured  for  steady-state  conditions.  The  flow  was  found  to 
be  supersonic  in  nature  with  an  imbedded  shock  for  all  but 
the  lowest  total  pressure  runs. 

The  analytic  method  used  one-dimensional  isentropic  and 
normal  shock  relations.  The  technique  provides  a  good 
"first-look"  at  the  nature  of  the  flow  and  is  compatible  for 
use  on  a  small  computer  systrem. 

The  experimental  data  was  obtained  by  using  an 
idealized  reed-valve  model.  The  data  acquistion  system  was 
computerized  and  used  a  scamvalve  system  to  measure  the 
pressures  m  the  valve  model.  Circular,  square,  and 
rectangular  valve  plate  geometries  were  investigated.  The 
inlet  was  circular  for  all  cases.  Flow- visualization  oil 
was  also  used  in  a  qualitative  approach  to  determine  the 
location  and  shape  of  the  shock  in  the  flow.  The  pressure 
measurements  and  the  flow  visualization  showed  the  shock  to 
be  circular  for  all  geometries  and  pressures  tested.  At 
high  total  pressure  runs  the  valve  plate  began  vibrating  and 
the  shock  location  became  smeared,  suggesting  that  the  shock 
was  also  oscillating. 


A  STUDY  UF  THE  FLUID-DYNAMIC  PRESSURE  FIELDS 
ON  COMPRESSOR  REED  VALVES 

I .  INTRODUCTION 

B  a  ckg r oun  d 

In  piston-type  compressors,  a  spring-metal  valve  plate 
is  often  used  lor  intake  and  exhaust  valves.  The  valve, 
hereafter  called  a  reed  valve,  has  wide  practical 
application  due  to  the  simplicity  of  operation  and 
manufacture,  and  the  minimal  number  of  moving  parts.  These 
qualities  normally  mean  long  life  and  as  such  these  valves 
are  in  common  usage  in  sealed  compressor  units  such  as 
refrigeration  systems.  These  systems  are  often  required  to 
operate  for  years  without  direct  maintenance. 

The  design  of  these  valves  has  historically  been  by 
trial  and  error.  Seemingly  small  design  changes  can  lean  to 
large  changes  in  the  operating  characteristics  of  the  valve. 
Rapid  and  violent  opening  and  closings  of  the  valve  during  a 
single  stroke  of  the  compressor  piston  after  slignt  design 
chanqes  to  previously  smoothly  operating  valves  have  been 
reported  il).  Modeling  of  the  valve  by  other  researchers  is 
often  from  a  structural  point  of  view.  In  other  work,  the 
pressure  on  the  face  of  the  valve  has  been  modeled  by 


directly  applying  the  total  pressure  of  the  flow  to  the  area 


ox  the  exposed  face.  The  present  thesis  work  investigated  a 
simple  X'luid-dynamic  model  of  the  pressure  field  for  valve¬ 
like  configurations.  The  pressure  field  on  the  up-flow  side 
of  the  reed  valve  was  of  primary  concern  in  this  study. 

The  thrust  of  this  work  was  to  develop  and  apply  an 
analytical  solution  to  valve-like  geometries  and  to  compare 
it  with  experimentally  acquired  data.  Certainly,  the 
problem  could  be  discretized  and  solved  using  an  elaborate 
computational  fluid  dynamics  approach,  but  the  intention 
here  was  to  take  a  "first  look”  with  the  simplest  approach. 
As  such,  the  analytical  solution  was  intended  to  be 
compatible  with  small  computer  systems. 

Initial  work  on  this  problem  was  accomplished  by 
Jumper  (2).  His  work  examined  the  valve  problem  from  a 
subsonic  compressible  flow  point  of  view.  One-dimensional 
isentropic  flow  relations  were  used  throughout.  While,  as 
will  be  shown,  the  flow  does  not  generally  remain  subsonic, 
insight  can  be  derived  from  his  findings.  One  such  insight 
is  the  fact  that  only  a  small  differential  between  the  inlet 
pressure  and  the  exit  pressure  is  required  to  choke  the  flow 
in  the  valve.  The  present  work  expanded  the  analytical 
treatment  of  J umper  to  include  supersonic  flow  and  the 
possibility  of  a  shock  wave  in  the  flow. 

Following  Jumper,  it  is  of  interest  to  recognize  the 
striking  geometric  similarity  between  the  reed-valve  problem 


*  *■  ■ 

AV. 


and  a  problem  involving  a  spool  of  thread  and  a  piece  of 


paper.  The  problem  involving  the  spool  of  thread  is  that  of 
a  child's  trick  where  a  pin  is  centered  in  a  piece  of  paper, 
then  the  paper  and  pin  are  placed  on  the  spool  of  thread. 
The  trick  is  to  blow  the  paper  off  of  the  spool.  The 
interesting  feature  of  the  trick  is  that  the  paper  cannot  be 
blown  off.  A  reed-valve  compressor  plate  and  a  reed-valve 
schematic  are  shown  in  Figures  1A  and  IB,  respectively.  The 
spool  trick  is  illustrated  in  Figures  2A  and  2B. 

The  geometric  similarity  between  the  spool  trick  and 
the  reed  valve  becomes  more  obvious  when  the  piece  of  paper 
is  trimmed  to  more  nearly  match  that  of  the  reed-valve 
configuration.  The  spool  trick  works  even  with  a  piece  of 
paper  which  barely  overlaps  the  hole  of  the  spool.  This 
simple  experimental  and  visual  aid  also  provided  a  starting 
point  for  the  design  of  the  experimental  apparatus. 

Other  flow  problems  with  similar  geometry  can  be  found. 
These  problems  are  typically  posed  as  incompressible  (3;  4) 
or  compressible  but  subsonic  as  in  the  case  of  the  radial 
diffuser  (5).  To  our  knowledge,  no  previous  work  has  been 
done  on  the  general  problem  which  was  the  focus  of  this 
study. 

Purpose  and  Scope 

The  purpose  and  scope  of  this  work  was  to  expand  on  the 
work  started  by  Jumper  and  to  compare  the  analytical 


solution  with  results  acquired  experimentally.  The  expanded 
analytical  solution  includes  subsonic  and  supersonic  theory 
with  the  possibility  of  normal  shocks  imbedded  in  the  flow. 
The  analytical  solution  of  the  problem  is  specifically 
taylored  to  be  preformed  on  a  small  computer  system.  To 
experimentally  validate  the  analytical  solution,  a  test  rig 
needed  to  be  designed  and  operated  over  a  range  of 
conditions  for  various  geometries  of  reed-valve  plates. 
Pressure  measurements  were  taken  for  the  various  geometries 


and  compared  with  the  analytical  solution. 


I I .  Theory 


Radial  Flow  Problem 


By  studying  the  geometry  of  the  spool  trick,  solution 
techniques  can  be  developed.  From  an  abstract  point  of 
view,  there  is  little  or  no  difference  between  the  physics 
of  the  spool  trick  and  the  reed-valve  problem  at  hand.  The 
two  problems  will  be  treated  here  as  the  same.  Figure  3 
details  the  problem  in  a  cross-sectional  view.  The  geometry 
of  the  problem  prescribes  a  flow  issuing  from  a  circular 
inlet.  The  flow  then  proceeds  radially  between  the  base 
plate  of  the  apparatus  and  a  free-floating  circular  valve 
plate.  Although  in  general  the  results  may  be  applied  to 
any  compressible  fluid,  it  will  be  assumed  that  air  is  the 
fluid  which  is  passing  through  the  valve. 

The  air  enters  this  idealized  reed  valve  from  a 
relatively  high  pressure  region.  This  high  pressure  region 
can  originate  from  the  action  of  a  piston,  as  in  the 
compressor,  a  high  pressure  source,  or  even  from  the 
pressure  exerted  by  the  operator  of  the  spool  trick  trying 
to  blow  the  piece  of  paper  off  of  the  spool.  When  the  air 
encounters  the  valve  plate,  it  is  turned  and  flows  radially 
to  the  exit  of  the  valve  plate.  In  many  situations  the  exit 
pressure  will  be  the  ambient  atmospheric  pressure  (though 
this  need  not  be  the  case).  Although  the  problem  at  hand 
will  be  treated  from  a  compressible  point  of  view,  a  look  at 


the  incompressible  analogue  is  helpful  for  a  background 


understanding  of  the  flow. 

Incompressible  Discussion 

Figure  4  shows  a  section  of  the  flow  which  is  helpful 
in  visualizing  the  problem  as  a  whole.  There  are  several 
cross-sectional  areas  of  the  flow  which  are  of  interest. 
The  first  is  the  area  of  the  inlet.  A  section  of  this  area 
is  identified  as  I  in  Figure  4.  The  next  sectional  area 
shown  is  that  of  the  annular  collar  around  the  inlet; 
identified  as  T.  Finally  the  exit  annulus  is  of  interest;  a 
section  of  which  is  shown  as  E.  The  areas  T  and  E  are 
dependent  on  the  height  <h)  of  the  valve  plate  above  the 
base  plate  of  the  apparatus.  A  radial  flow  of  this  type, 
for  an  incompressible  fluid,  is  similar  to  a  two-dimensional 
potential  source  flow  (4:134).  Now,  with  a  rough  idea  of 
the  geometry  through  which  the  air  flows,  equations  for 
incompressible  f low  through  a  variable- area  stream  tube  may 
be  applied. 

In  the  incompressible  case,  pressure  and  velocity  are 
related  by  the  incompressible  form  of  Bernoulli's  equation 

Po  =  P  ‘  (1/2)  p  V'  (  1  > 

where  P  is  the  pressure,  p  is  the  density,  and  V  is  the 
velocity.  The  subscript  "0 "  refers  to  total  or  stagnation 


According  to  continuity,  the  mass  flow  rate 

m  =  pAV  =  constant  (2) 

can  be  rearranged  and  differentiated  to  yield  a  relation  for 
velocity  and  area  in  an  incompressible,  steady  flow,  so  that 

1  dV  =  -  1  dA 

V  dx  A  dx  ( 3 ) 

This  equation  demonstrates  that  for  subsonic,  incompressible 
flow,  an  increase  in  the  area  of  a  given  stream  tube  is 
accompanied  with  a  decrease  in  the  velocity  of  the  fluid  in 
the  stream  tube.  Similarly,  a  decrease  in  area  yields  an 
increase  in  velocity.  For  the  problem  of  the  radial  flow 
from  the  inlet  annulus  to  the  exit  annulus,  the  area  is 
increasing  which  indicates  that  the  fluid  velocity  must 
decrease  in  the  radial  direction. 

The  radial  decrease  in  velocity  can  also  be  identified 
with  an  increase  in  the  static  pressure.  Referring  back  to 
Eq  (1),  a  decrease  in  velocity  for  an  incompressible  fluid 
must  be  accompanied  with  an  increase  in  static  pressure  for 
a  fixed  total  pressure.  Finally,  for  subsonic  flow,  the 
exit  pressure  must  equal  the  ambient  pressure  at  the  exit, 
this  being  due  to  the  ability  of  isentropic  pressure  signals 
to  comunicate  upstream  at  the  speed  of  sound.  This 

communication  is  unavoidable  because  the  flow  is  subsonic. 
The  ambient  pressure  signal  at  the  exit  is  able  to 


communicate  upstream  into  the  valve  and  modify  the  flow 
until  the  exit  pressure  exactly  matches  the  pressure  of  the 
ambient  environment. 

The  height  of  the  valve  plate  above  the  base  plate,  h, 
is  an  interesting  parameter  of  the  valve  problem.  Because  h 
is  a  factor  in  determining  the  area  of  the  annular  inlet 
collar  (see  Figure  4)  it  can  dramatically  alter  the  nature 
of  the  flow.  If  h  is  small,  then  the  area  of  the  inlet 
collar  can  be  less  than  that  of  the  iniet.  This  would 
effectively  produce  a  channel  area  distribution  which  is 
similar  to  a  convergent -divergent  nozzle.  A  bit  of  reverse 
logic  shows  that  indeed  this  is  the  only  reasonable  physical 
case  to  expect  for  a  free-floating  valve  plate. 

For  a  moment,  let's  assume  that  h  is  sufficiently  large 
to  insure  that  the  inlet  collar  area  (cf.  Figure  4  and 
sectional  area  T)  is  greater  than  the  area  of  the  inlet. 
Also,  assume  that  the  flow  is  not  allowed  to  separate.  For 
How  to  be  in  the  desired  direction,  the  inlet  (reservoir) 
pressure  must  be  greater  than  the  ambient  pressure, 
P«ti  >  P,  .  As  the  flow  passes  the  inlet,  the  area 
increases  and  continues  to  increase  to  the  exit.  As  shown 
before,  this  increase  in  the  area  of  the  flow  channel  is 
accompanied  with  a  decrease  in  velocity  and  an  increase  in 
static  pressure.  This  hypothetical  case  leads  to  the  exit 
pressure  being  greater  than  the  reservoir  pressure  or 
Pt  >  P« i •  >  P«  Recall  however,  that  for  a  subsonic  flow. 


the  exit  pressure  must  equal  the  ambient  pressure. 
There! ore,  the  exit  pressure  cannot  be  greater  than  that  oi 


the  ambient  pressure  as  described  in  this  hypothetical  case. 
This  descrepency  shows  that  a  large  h  is  not  compatible  with 


a  free  floating  valve  plate.  Since  the  flow  will  simply 
adjust  until  P*  =  P„  ,  the  proper  sequence  of  pressures 
must  be  Pm  >  Pe  =  P„  .  This  sequence  occurs  when  h  is 
small  enough  to  make  the  inlet  collar  area  less  than  the 
inlet  area;  the  inlet  collar,  in  effect,  throttles  the  flow 
in  the  valve. 

The  height  of  the  valve  plate  is  a  function  of  the  mass 


flow  rate,  as  defined  by  continuity  and  the  geometry  of  the 


inlet  collar. 


The  area  of  the  collar  is  simply 


Ac  =  2nR,<Ltrh  .  Combining  this  with  Eq  (2>  completes  a 
full  specif ication  of  the  incompressible  problem. 

This  simplif led  incompressible  look  at  the  problem 
shows  why  the  operator  of  the  spool  trick  is  unable  to  blow 
the  paper  off  of  the  spool.  The  entire  radial  flow  between 
the  inlet  collar  and  the  exit  collar  is  at  a  pressure  below 
that  of  the  ambient  environment.  The  entire  external  side 
of  the  paper  is  acted  upon  by  ambient  pressure.  This  force 
imbalance  is  opposed  by  the  force  necessary  to  effect  a 
momentum  change  in  the  flow  to  turn  the  corner  at  location  T 
t see  Figure  4;  and  possibly  the  weight  of  the  paper.  For 
geometries  of  the  trick  where  the  piece  of  paper  is 
approximately  the  same  size  as  the  spool  of  thread,  the 
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operator  is  physically  unable  to  produce  enough  pressure  to 
dislodge  the  paper.  It  will  be  shown  that  the  flow  enters 
the  compressible  regime  and  chokes  at  location  T  far  before 
the  force  due  to  the  momentum  change  can  dislodge  the  paper. 


Compressible  Solution 

Many  of  the  basic  features  of  the  flow  as  identified  in 
the  incompressible  discussion  still  hold.  The  concept  of 
the  area  of  the  channel  being  similar  to  a  convergent- 
divergent  nozzle  is  valid,  but  the  mathematical  decnption 
in  Eq  (3)  will  require  modification  for  compressible  flow. 
With  a  proper  pressure  differential  between  the  reservoir 
and  the  ambient  environment,  the  flow  may  choke  at  the  inlet 
collar  and  it  is  possible  to  have  a  supersonic  section  of 
the  flow.  A  look  at  one-dimensional  isentropic  flow 
conditions  shows  why  an  increasing  area  accelerates  the  flow 
supersonically.  The  differential  form  of  the  continuity  and 
momentum  equations  for  one-dimensional  flow  is  given  by 


VdV  »  dP  =0 

P 


If  we  assume  isentropic  flow  and  incorporate  the 
definition  of  the  speed  of  sound,  ( dP/ds ; %  ) *  ■  3  ,  and  the 
continuity  equation,  Eq  (4)  may  be  manipulated  <b:209-213) 
to  yield 


(  Ma  -  1  )  1  dV  =  1  dA 
V  dx  A  dx 


Three  cases  must  be  considered. 

For  M  <  1  ,  dV/dx  and  dA/dx  are  opposite  in  sign. 
This  is  similar  to  Eq  (3).  Hence  a  subsonic  flow 
decelerates  in  a  diverging  channel  and  accelerates  in  a 
converging  channel. 

For  M  >  1  ,  dV/dx  and  dA/dx  have  the  same  sign. 
Hence  a  supersonic  flow  accelerates  in  an  expanding  channel 
and  decelerates  in  a  converging  channel.  Thus  supersonic 
behavior  in  a  channel  is  opposite  to  that  of  subsonic 
behavior. 

For  II  =  1  ,  the  most  appropriate  consideration  from 
physical  grounds  is  that  dA/dx  =  0  .  We  can  conclude  that 
the  sonic  condition  can  only  occur  at  a  constriction  or 
throat. 

Although,  it  may  at  first  appear  that  the  problem  at 
hand  is  not  suited  to  a  one-dimensional  analysis,  the 
relative  simplicity  of  the  one- dimensional  i sen tropic 
equations  makes  them  attractive.  Further,  their  use  is 
compatible  with  a  basic  tenent  of  this  work,  to  use 
techniques  which  are  compatible  with  small  computer  systems. 
The  reed-valve  problem  presents  a  geometry  to  the  flow  which 
is  comparable  to  a  convergent -divergent  nozzle.  The  one- 
dimensional  lsentropic  relations  are  applicable  to  a  wide 
variety  of  nozzle-type  problems,  with  reasonable  results, 
even  when  the  one-dimensionality  of  the  problem  is 
questionable.  This  work  will  use  these  relations  and  the 


normal  shock  relations  for  the  analytical  solution  to  the 
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problem,  and,  as  will  be  shown,  this  method  yields 
reasonably  good  agreement  with  experiment. 

Numerical  Solution 

A  computer  program  was  written  to  solve  the  reed-valve 
problem.  The  program,  SOLVE,  is  listed  in  Appendix  A. 
SOLVE  was  written  in  BASIC  and  uses  the  assumption  that  the 
flow  is  steady,  one-dimensional,  at  a  minimum  piece-wise 
isentropic  and  has  the  possibility  of  imbedded  snocks.  The 
flow  is  assumed  to  be  air  and  was  treated  as  a  perfect  gas. 
The  program  was  run  on  a  Columbia  Data  Products  personal 
computer  (IBM  compatible)  which  has  256K  RAM.  SOLVE 

produces  a  solution  to  the  reed-valve  problem  in 
approximately  60  seconds. 

The  anaylitical  solution  includes  a  number  of 
isentropic  flow  parameters  (7;  8:47-75).  In  the  following 
discussion  the  "0 "  subscript  again  will  refer  to  total  or 
stagnation  conditions  for  the  specific  parameter.  Equations 
tor  the  parameters  of  interest  have  been  reduced  to  forms 
which  require  input  of  the  total  condition  of  the  parameter 
and  the  local  Mach  number.  A  function  oi  the  local  Mach 
number  is  a  common  factor  to  many  of  these  calculations  and 
is  given  by 

=  C  1  ♦  (  7  - 1  )  M4  /  2  ] 
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This  function  of  Mach  number  was  then  used  to  determine 
the  following  properties: 

Pressure  was  obtained  by 


P  ■  _ Efl _ 

r  7/0-1) 

|  FNBOB(M)  | 
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Density  was  found  by  applying 

P  =  _ Pji _ 

r  i  1/(  -*> 

j  FNBQB(M)  ; 

i-  -i  (8) 

The  local  temperature  was  found  using 

T  =  _ Tq _ 

FNBOB <  M )  (9) 


Also,  using  the  total  temperature,  it  was  convenient  to 
define  a  total  sonic  velocity 
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So  that  the  local  sonic  veloci t y  could  be  found  using 


c  =  _ ce _ _ _ 

~  r .  ”  ~  i  1/2 

j  FNBOB ( M  >  j 


ill) 


The  local  velocity  of  the  flow  was  now  found  by 


In  order  to  find  the  Mach  number  at  each  location  of 


interest  two  equations  relating  the  Mach  number  to  the 
pressures  and  Mach  number  to  area  ratios  were  required. 
First,  for  an  isentropic  flow,  the  Mach  number  at  the  exit 
could  be  found  by  knowing  the  pressure  ratio  between  the 
total  pressure  and  the  exit  pressure.  The  relation  is 


M*  = 
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Once  the  exit  Mach  number  was  known,  the  Mach  number  at 
any  other  location  could  be  determined  from 
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It  is  important  to  note  that  while  the  radius  ratio  is  used 
in  the  preceedmg  equation;  the  usual  ratio  is  that  of  the 
area  ratios.  For  the  circular  reed-valve  model,  the  radii 
ratio  is  identically  the  area  ratio  from  the  inlet  collar  to 
the  exit  collar.  Recall  that  the  area  of  the  inlet  collar 
is  Ac  =  2rfii »L  t i h  and  the  area  of  the  exit  collar 
is  Ac  =  2nRt  h  .  The  common  factor  of  2nh  cancels  from 
the  ratio  of  the  two  areas.  It  is  emphasized  here  that  this 
problem  is,  therefore,  fully  specified  from  this  theoretical 
point  of  view  without  the  variable  h  entering  the  problem. 


If  a  shock  occurs  in  the  flow 


several  parameters 


undergo  a  change  which  was  modeled  by  a  jump  discontinuity . 
The  circular  geometry  of  the  reed  valve  model  and  the  radial 
nature  of  the  flow  through  the  valve  suggest  that  any  shock 
in  the  flow  will  be  normal  to  the  flow.  Although  the  shock 
will  be  circular  when  viewed  from  above,  the  flow  sees  a 
simple  normal  shock,  thus  the  normal  shock  relations  are 
appropriate  for  this  case. 

The  subscripts  "  I  H  and  "»  "  denote  the  upstream  and 
downstream  values,  respectively,  of  a  parameter  across  the 
shock.  If  the  upstream  Mach  number  is  known,  the  downstream 
Mach  number  is  given  by 


The  total  pressure  decreases  across  the  normal  shock. 
Total  pressure  is  the  parameter  which  is  used  in  SOLVE  to 
track  the  changes  in  the  other  parameters  after  the  shock . 
Static  pressure  in  particular  is  of  interest  because  the 
Ilow  after  the  normal  shock  will  be  subsonic  and  therefore 
the  static  pressure  at  the  exit  must  match  the  ambient 
pressure.  This  was  the  parameter  to  be  matched  when 
determining  the  location  of  the  shock.  The  total  pressure 
on  the  downstream  side  of  the  shock  is  given  by 
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With  the  background  of  the  specific  equations  in  place, 
a  dicussion  of  the  logic  of  SOLVE  follows.  An  overview  of 
the  entire  program  is  given,  followed  by  additional 
overviews  of  specifics  of  the  subsonic  and  the  supersonic 
solution  techniques. 

Figure  5  is  a  flow  chart  of  SOLVE.  Unless  otherwise 
specified  the  program  defaults  to  standard  day  conditions 
for  temperature  and  pressure.  The  inlet  radius  for  the 
valve  is  also  set  to  default  to  the  value  of  the  actual 
reed-valve  model  tested.  An  alternate  input  section  is 
interacti vely  available  to  the  user  to  allow  specific  inputs 
to  match  a  specific  problem.  The  default  option  is  useful 
tor  a  general  understanding  of  a  preliminary  problem; 
however,  it  is  important  to  use  the  alternate  input  section 
it  one  wishes  to  compare  theoretical  values  with 
experimental  data.  Next,  following  the  input  section,  SOLVE 
sets  the  definitions  of  various  parameters,  descretizes  the 
radial  domain  of  the  flow  field,  and  determines  the  nature 
of  the  problem. 

Depending  on  the  ratio  of  reservoir  pressure  to  the 
exit  (ambient)  pressure  the  nature  of  the  problem  can  be 
subsonic,  supersonic  with  imbedded  shocks,  or  fully 
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supersonic.  The  main  distinction  which  SOLVE  is  concerned 
with  is  whether  the  problem  will  require  a  subsonic  or  a 
supersonic  solution  technique.  The  program  uses  Eq  (14)  and 
a  fixed  point  iteration  scheme  (9:88-95)  to  find  the  Mach 
number  at  the  exit  for  the  case  where  the  flow  is  exactly 
sonic  at  the  throat  (inlet  collar)  but  subsonic  elsewhere. 
This  condition  will  be  refered  to  as  the  star  condition  for 
the  problem.  The  flow  is  isentropic  throughout  the  field 
tor  this  case.  Knowing  the  exit  Mach  number  for  the  star 
condition  the  pressure  ratio  for  the  star  condition  can  be 
found  using  Eq  (7).  By  comparing  the  pressure  ratio  of  the 
problem  with  the  star  condition  pressure  ratio,  the  nature 
of  the  problem  is  established.  Following  the  appropriate 
solution  technique,  SOLVE  outputs  the  data  to  the  screen  for 
viewing  and  will  save  pressure  or  pressure  ratio  data  to  a 
disk  file.  The  files  are  formatted  for  easy  use  with 
graphics  software. 

It  the  specif led  pressure  ratio  is  less  than  the 
pressure  ratio  for  the  star  condition,  the  problem  will  be 
entirely  subsonic.  A  flow  chart  of  the  subsonic  solution  is 
presented  m  Figure  6.  Here  the  exit  Mach  number  is  found 
using  Eq  (  L  3 ) .  The  Mach  numbers  at  the  incremental 
locations  are  then  found  using  Eq  (14)  and  fixed  point 
iteration.  Finally,  the  remainder  of  the  properties  at  the 
incremental  locations  are  found  using  Eqs  (7-12).  The 


program  branches  to  the  output  section  with  the  subsonic 
description  of  the  flow. 

If  the  flow  is  supersonic  in  nature,  the  program 
assumes  that  the  flow  also  has  an  imbedded  normal  shock. 
Schreier  (8:75)  gives  a  trial-and-error  method  for  solving 
this  problem.  Some  of  his  procedure  has  been  modified  to 
make  it  more  applicable  to  a  computer  program  and  the  reed- 
valve  problem.  A  flow  chart  of  the  supersonic  solution  is 
shown  in  Figure  7.  The  steps  used  in  SOLVE  are: 

1.  Postulate  an  upstream  Mach  number  for  fixed  shock. 

2.  Calculate  the  radius  ratio  using  Eq  (14). 

3.  Find  the  downstream  Mach  number  using  Eq  (15). 

4.  Find  the  downstream  total  pressure  using  Eq  (16). 

5.  Find  the  exit  Mach  number  using  Eq  (13). 

6.  Find  the  exit  pressure  using  Eq  (7). 

7 .  Compare  this  pressure  with  the  ambient  c:  back 

pressure.  If  they  are  equal  to  within  a  prescribed 


tolerance 

the 

shock  is 

in  the 

desired  location. 

Otherwise, 

the 

process  is 

repeated 

by  changing  the 

position  of  the  shock.  This  is  done  by  postulating  a 
new  upstream  Mach  number. 

A  marching  technique  is  used  to  place  the  shock  in  the 
proper  location.  After  the  initial  location  is  evaluated, 
the  program  takes  a  step  increase  or  decrease  in  the 
upstream  Mach  number  and  re-evaluates  the  problem.  This  is 
repeated,  stepping  in  the  same  direction,  until  it  has 
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overstepped  the  solution. 


At  that  time  the  march  changes 


direction  and  cuts  the  size  of  the  step  in  half.  The  march 
continues  m  the  new  direction  until  the  solution  is  again 
overstepped.  This  whole  process  continues  until  the  exit 
pressure  matches  the  prescribed  back  pressure  to  within  a 
specified  tolerance.  This  marching  technique  is  also  used 
to  determine  the  Mach  number  at  specified  area  ratios  for 
the  supersonic  portion  of  the  flow  field  because  the  fixed 
point  iteration  scheme  does  not  converge  to  the  supersonic 
value  of  Mach  number.  Once  the  shock  has  been  determined  to 
be  in  the  proper  location,  the  program  delivers  the 
incremental  properties  to  the  output  section. 

A  collection  of  theorectical  pressure  ratio  versus  the 
radial  location  data  from  SOLVE  is  plotted  in  Figure  8. 
Figure  8  includes  a  subsonic  curve,  a  nearly  star  condition 
curve,  and  several  supersonic  curves  with  imbedded  shocks. 
One  feature  to  keep  in  mind  about  Figure  8  is  that  the 
l sen tropic  preshock  port ions  of  the  supersonic  curves  lay  on 
top  of  each  other  up  to  the  point  where  they  shock.  This 
will  be  a  pattern  to  be  looking  for  wnen  viewing  the 
experimental  data.  The  theoretical  curves  presented  here 
are  similar  to  those  in  texts  describing  the  liow  in  a  Laval 
nozzle  i 6 : 21 1 -213, 242 ;  10:127-130)  because  the  geometry  of 
the  flow  in  the  reed  valve  is  similar  to  a  convergent- 
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divergent  nozzle.  It  should  be  noted,  however,  that  the 
convergence  in  the  reed-valve  problem  is  a  very  sharp 
decrease,  basically  a  step,  and  the  divergent  section  has  a 
rapidly  increasing  area.  An  equivalent  two-dimensional 
nozzle  would  have  a  divergence  half  angle,  A,  of 
approximately  72.3  degrees  (  tan  A  =  n  ).  For  a  circular 
duct  nozzle,  the  applicable  divergence  half  angle 
relationship  is  tan  A  =  (2h/Ri  )*•  5  .  It  is  emphasized 
here  that  the  reed-valve  geometry  does  not  match  either  the 
two-dimensional  or  three-dimensional  equivalent  nozzle,  but 
a  comparison  with  these  "standard"  nozzles  is  interesting. 

Other  Theoretical  Considerations 

It  is  interesting  to  look  at  what  range  of  pressure 
ratios  that  cause  the  reed-valve  or  the  spool  problem  to 
change  from  a  subsonic  nature  to  a  supersonic  problem. 
Assume  an  inlet,  circular  in  cross  section,  with  a  radius  of 
1/6-inch  and  a  reed-valve  plate  (or  piece  of  paper)  with  a 
radius  of  1-inch  (representative  of  the  spool  trick).  Hy 
using  SOLVE,  it  can  be  seen  that  the  problem  changes 
character  at  a  pressure  approximately  0.  0b5  psi  above 
ambient  (standard  day  assumed).  This  seems  incredibly  low 
and  means  that  lor  inlet  pressures  greater  than  this,  the 
flow  will  have  a  supersonic  section  with  an  imbedded  shock. 
When  performing  the  spool  trick,  an  average  adult  can 
maintain  a  pressure  on  the  order  of  1  to  2  psi.  Obviously, 


even  the  spool  trick  must  be  treated  as  a  compressible 
supersonic  problem. 

A  more  formal  presentation  of  what  pressure  ratios  are 
required  to  theoretically  produce  a  supersonic  problem  is 
given  in  Figure  9.  An  example  here  is  also  instructive. 
For  a  radius  ratio  of  0. 8  and  the  exit  pressure  at  ambient, 
the  valve  would  have  sonic  flow  at  the  throat  for  a  gauge 
pressure  of  approximately  3. 6-psig.  This  example 
corresponds  to  a  valve  plate  with  a  radius  of  0. 344-inches 
for  an  inlet  of  radius  of  0. 281 - inches.  The  overlapping  lip 
of  the  valve  plate  over  the  inlet  hole  then  is  0.063-inches. 
This  amount  of  overlapping  for  a  valve  lip,  about  1/16  inch, 
is  representative  of  the  order  of  the  measurements  for  an 
actual  reed-valve  configuration  (see  Figure  1). 

Although  the  theoretical  treatment  of  the  reed-valve 
problem  is  inviscid,  boundary  layers  should  be  expected  to 
be  present  in  the  valve.  When  the  air  initially  impinges  on 
the  valve  plate  a  stagnation  or  Hiemenz  How  (11:95-99)  is 
likely.  Further  into  the  channel  boundary  layer  growth  for 
a  divergent  channel  could  be  expected  (11:106,223).  The 
boundary  layer  will  have  a  displacement  thickness  which 
should  have  an  effect  on  the  experimental  results.  The 
height  of  the  valve  plate  above  the  base  plate,  if  constant, 
is  not  a  required  parameter  in  the  theoretical  solution  of 
the  reed-valve  problem;  however,  boundary-layer  growth  could 


Radius  Ratio  vs.  Pressure  Ratio  From 
Analytical  Solution 


change  the  etiective  channel  height  from  a  constant  to  a 
variable.  For  an  increasing  displacement  thickness,  the 
area  of  the  channel  would  not  grow  as  rapidly  in  the  viscous 
case  as  in  the  theoretical  case.  One  might  assume,  then, 
that  the  pressure  ratios  at  a  specific  location  would  have  a 
value  corresponding  to  a  smaller  area  ratio.  One  would  also 
expect  other  considerations,  such  as  separation  of  the  flow, 
to  effect  the  solution. 
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gure  11.  Schematic  of  Work  Apparatus 
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.'ne  air  wmcn  t*awed  to  tne  test  nq  was  metered 
tnrouqn  a  regulator.  Two  air  sources  were  avai*aoie  tor  tne 


n 


testing.  The  tirst  was  the  tmnaing  air  supply  ana  tne 
second  was  bottled  air.  The  building  air  was  used  for  data 
runs  up  to  about  Big  psig,  bottlea  air  was  usea  for  Higher 
pressures.  The  total  pressure  of  the  reservoir  which 
suppuea  tne  air  was  measured  on  eitner  a  iiln-incn  mercury 
•nanometer  or  a  irild-incn  oouraon-type  pressure  gauge.  A 
complete  listing  oi  tne  equipment  uses  in  tne 
experimentation  is  presented  in  Appendix  B. 


Test  Kig 

The  test  rig  used  for  the  experimental  investigation  is 
shown  in  Kigure  IsA  and  lib.  As  stated  ear.ier,  the  design 


of  tne  test  rig  ioiiowea  Irom  considering  tr.e  r  eea  vai  ve 
,-igure  11  ana  tne  spool  tncx  (Figure  i  >  .  The  similarity 
ci  the  spooi  trie;-:  ana  tne  test  rig  is  again  easny  seen. 
iacn  provide  a  channel  lor  tne  now  wmcn  expands  radiai.y 
outwarc.  ;  he  channel  is  constrainea  oy  a  noating  valve 
ciate  above  ana  tne  rixea  oase  piate  oeiow,  uue  to  tne  test 
.  .  c  teing  more  adaptao.e  to  being  connected  to  a  nigh 


pressure  reservoir,  no  lormai  measurements  were  taxen  using 


ipuo i  ui  t  n i eaa. 


ne  test  rig  was  lacncated  oy  Aril'  moaei  snap 
personnel  irom  aiuminum  stoex.  ine  w*4-incn  diameter  1  me t 
nc.e  was  centered  m  the  oase  plate.  Fortv  surtace  pressure 


V  «* 


PRESSURE  (INCHES  Hg) 
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Test  Procedure 

A  computer  program,  TEST,  was  written  to  direct  the 
acquisition  of  the  pressure  measurements  in  the  reed-valve 
model.  The  actual  test  procedure  was  straight  forward.  The 
equipment  was  turned  on  and  allowed  to  warm  up  and 
stabilized  for  approximately  15  minutes.  During  the  warm-up 
period  barometric  pressure  and  ambient  temperature  readings 
were  made.  After  warming  up,  the  bridge  excitation  voltage, 
which  was  sent  from  the  power  supply  to  the  scamvalve,  was 
checked  and  adjusted  as  necessary.  A  calibration  run  for 
the  transducer  was  accomplished  next.  This  calibration  was 
done  lor  each  set  of  tests,  but  it  is  noted  here  that  there 
was  no  measurable  difference  in  the  calibration  curve  over 
the  entire  duration  of  the  testing.  Once  the  calibration 
was  checked,  the  testing  was  ready  to  procede. 

The  configuration  to  be  tested  was  selected  and  set  up 
on  the  test  rig.  The  inlet  pressure  was  set  by  adjusting 
the  regulator  on  the  inlet  line.  The  system  was  again 
allowed  to  stabilize.  While  stabilizing,  the  program  TEST 
was  started  with  a  run  command.  TEST  interactively 
requested  the  test  condition  and  then  awaited  a  command  to 
commence  gathering  data.  When  the  system  haa  stabilized  the 
operator  input  a  line  ieed  (return)  command  and  PEbT 
directed  the  gathering  of  the  data.  The  data  runs  lasted 
approximately  30  seconds.  During  the  data  run,  ail  40 


K 

V  ■ 


suriace  pressure 


ports  were  sampled. 


Each  port  was 


individually  sampled  20  times  before  the  scam valve  was 
stepped  to  the  next  port.  The  value  of  the  reading  for  each 
port  measurement  was  the  average  of  the  multiple  sampling. 
Upon  completion  of  the  data  run,  TEST  could  be  made  to  save 
the  data  to  a  disk  file. 

This  procedure  of  testing  continued  until  the  operator 
was  finished  with  the  desired  runs  for  the  day.  Finally, 
the  barometric  pressure  and  temperature  were  again  measured. 

F low  Visualiz at i on 

In  order  to  obtain  a  more  specific  idea  of  the  shock 
location  and  shape,  a  flow  visualization  technique  was 
devised.  Wind  tunnel  flaw  visualization  oil  was  used.  Due 
to  the  oil  mixture  containing  white  pigment,  the  test  rig 
was  painted  black  to  insure  better  contrast  for  photography. 
The  oil  was  dotted  on  the  base  plate  of  the  test  rig  as 
shown  in  Figure  16A.  In  order  to  keep  the  oil  undisturbed 
(except  by  the  flow)  the  tests  were  run  with  the  rig  up-side 
down  as  shown  in  Figure  lbB.  This  required  a  slightly 
different  configuration  for  the  test  rig.  The  legs  were 
reversed  in  their  holes  and  knife  edge  supports  replaced  the 
centering  pins  (cl.  Figures  12  and  lbC;.  After  the  oil  was 
in  place  tne  rig  was  inverted  and  the  air  was  turned  on  at  a 
moderate  rate.  The  valve  plate  was  then  raised  slowly 
toward  the  inlet.  Recall,  that  the  problem  is  sell 
adjusting  for  the  height  of  the  valve  plate,  so  once  the 
plate  is  close,  the  physics  of  the  problem  takes  over  and 
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snatches  the  plate.  Alter  the  valve  plate  is  captured  by 
the  test  rig,  the  reservoir  pressure  was  adjusted  to  the 
test  condition.  The  test  pressure  was  maintained  lor 
approximately  20-seconds.  The  air  supply  was  then  cut  off 
and  the  valve  plate  fell  to  the  table,  leaving  the  oil  flow 
undisturbed.  The  rig  was  then  turned  back  over  for 
photographs  of  the  flow  pattern. 


V. 


Data  Reduction 

l'he  data  were  saved  to  a  disk  iiie  during  the  test 
runs.  These  data  files  were  accessed  by  a  data  reduction 
program,  REDDATA.  REDDATA  was  written  in  BASIC  and  is 
listed  in  Appendix  D.  The  program  took  the  raw  data  which 
was  collected  by  TEST,  applied  the  appropriate  corrections 
to  the  data,  converted  the  pressure  readings  to  absolute 
pressure  readings  in  psi,  and  then  saved  the  reduced  data  to 
fij.es  for  plotting. 

The  data  which  was  collected  by  TEST  was  presented  in 
millivolts.  REDDATA  took  the  millivolt  readings  and  used 
the  calibration  curve  to  convert  the  millivolt  readings  to 
the  equivalent  manometer  reading  in  inches  of  mercury.  The 
readings  required  correction  at  this  point.  The  equivalent 
manometer  reading  was  corrected  for  temperature  and  gravity, 
l’he  temperature  correction  accounted  for  thermal  expansion. 
The  gravity  correction  accounted  for  the  elevation  and 
latitude  of  the  test  site.  Tables  for  the  corrections  are 
available  in  several  experimental  guides  or  meteorological 
texts  (ID;  14:13,  13B-2<i2  )  . 

l'he  pressure  readings  to  this  point  were  in  inches  of 
mercury.  To  get  an  absolute  pressure,  the  ambient  pressure 
had  to  be  added  to  the  data.  The  ambient  barometric  reading 
was  taken  at  test  time  using  a  Fortin-type  barometer.  This 
reading  was  also  in  inches  of  mercury  and  had  to  be 
corrected.  Once  corrected,  the  pressure  readings  could  be 


added  and  converted  to  psx  readings. 


Similanly,  the  total 


(reservoir)  pressure  reading  was  corrected  and  converted. 

REDDATA  took  the  absolute  pressure  readings  and 
presented  the  data  in  a  format  similar  to  that  which  SOLVE 
used  for  the  theoretical  values.  REDDATA  used  the  symmetry 
of  the  problem  to  give  an  overall  view  of  the  flow  field 
properties  in  the  reed-valve  test  rig.  For  the  circular 
valve  plate  the  flow  field  along  each  radial  arm  was 
equivalent.  REDDATA  averaged  the  eight  radial  profiles  to 
get  the  average  description  of  the  radial  flow -field 
pressure  under  the  valve  plate.  The  reducted  data  was  then 
saved  to  a  disk  file  for  plotting.  REDDATA  could  save  the 
data  in  a  pressure  ratio  format  or  in  an  absolute  pressure 
xormat. 

Variations  of  RED DAT  A  were  used  to  reduce  the  data  tor 
the  square  and  rectangular  valve  plates.  The  square  plate 
nad  two  sets  01  symmetric  profiles.  The  radial  arms  which 
were  parallel  to  the  V)  and  9M-degree  axes  of  the  square  were 
symmetric  and  could  oe  averaged.  The  *4b  and  -4b-degree 
diagonals  were  also  symmetric.  For  the  rectangle,  3  sets  of 
symmetric  profiles  are  used.  Figure  17  identifies  tne 
radial  arms  which  display  symmetry  lor  the  various  plates. 
The  variations  to  the  REDDATA  program  required  to  present 
the  data  for  the  square  and  rectangular  disks  are  straight 
forward  and  not  presented. 


A  nnal  program  was 


written  to 


compare 


the  area 


presen tea  by  the  reed- valve  test  rig  with  the  equivalent  or 
effective  area  of  an  ideal  nozzle.  This  program,  EKFAREA, 
is  a  modification  of  SOLVE.  The  purpose  of  writting  EFFAREA 
was  to  look  at  various  area  ratio  profiles  in  order  to  find 
a  best-fit  match  of  the  data.  The  modifications  to  SOLVE 
consisted  of  replacing  all  references  to  radii  ratios  with 
area  ratios  and  by  specifying  specific  area  ratios  at  all 
incremental  locations.  Two  methods  were  used  to  define  the 
ef fective  area  ratios. 

The  first  profile,  derived  by  method  1,  was  formed  by 
fairing  in  the  curve  for  the  data  which  was  analogous  to  the 
supersonic  portion  of  the  theoretical  curve.  Using  the 
taired  curve,  pressures  ratios  could  be  projected  back  to 
tne  theoretical  curve  and  effective  area  ratios  could  be 
determined. 

The  second  profile,  using  method  2,  was  basically 
similar  to  the  profile  used  in  SOLVE.  However,  the  area  of 
the  throat  was  varied  while  leaving  the  other  incremental 
areas  unchanged.  This  scheme  may  be  similar  to  a  profile 
wmcn  was  due  to  a  blockage  at  of  near  the  corner  which 
defined  tne  throat.  The  location  of  the  shock  was  used  as 
tne  comparison  criterion  to  determine  what  factor  should 
multiply  the  throat  area.  After  comparing  several  pressure 
profiles,  the  factor  which  multiplied  the  throat  are  was 
determined  to  be  0.64. 


VI. 


Kesuiis  ang  uiscussion 


uenerai 

Data  from  7<D  test  runs  were  saved  to  disk  files  for 
reduction.  Tests  on  three  configurations  were  accomplished . 
In  this  section  the  configurations  will  be  referred  to  as 
configuration  1,  2,  and  3.  The  first  coni igur at ion  was  the 
circular  plate,  the  second  was  the  square  plate,  and  the 
third  was  the  rectangular  plate.  Configuration  1  was  tested 
in  33  runs;  configuration  2,  31  runs;  and  configuration  j, 
o  runs.  Testing  was  accomplished  an  S  separate  days. 
Duplicate  conditions  were  tested  on  different  days  to  check 
the  repeatability  of  the  test  data.  The  repeatability  for 
all  checks  was  excellent.  The  pressure  ratio  data  points 


for  similar  runs  from  different  days  were  basically 


coincident. 


Configuration  I 

Test  runs  lor  the  circular,  2-inch  diameter  valve  plate 
were  performed  for  gauge  pressures  from  1.2-inches  of 
mercury  to  21VD-inches  of  mercury.  Figure  la  shows  a 
re presen  tative  sampling  of  5  data  runs.  A  point  to  note 
from  tne  figure  is  that  port  ions  of  tnese  curves  closer  y 
overlay  one  another.  i'his  is  very  representative  of  the 
lsentropic  subsonic  to  supersonic  branch  of  the  one- 
aimensional  theoretical  analysis  (cl.  Figure  U ) .  The  exact 


locations  of  the  normal  shocks  m  these  ilows  are  unknown. 


PRESSURE  RATI 


0.00  0.20  0.40  0.60  0.30  1.00  1.20 

RADIAL  DISTANCE  (INCHES) 

□  15. 88  PSIA  A  20.46  PSIA  0  28.96  PSIA  X  72.90  PSIA  +104.9  PSIA 


Figure  18.  Representative  Data;  Pressure  Ratio  vs. 

Radial  Location;  Five  Runs 


it  is,  however,  iairiy  clear  Iron  the  jump  in  the  pressure 
ratio,  that  a  shock  must  lay  between  certain  data  points. 
Data  points  on  opposite  sides  of  the  shock  are  connected  by 
a  dotted  line.  Keep  in  mind  that  the  shock  caused  a  steeper 
jump  in  pressure  ratio  than  this  dotted  line  represents. 

Figures  19-28  are  direct  comparisons  of  analytical  and 
experimental  data.  Figure  29  is  a  compilation  of  several 
pressure  ratio  plots.  The  theorect ical  curve  for  a  specif ic 
reservoir  to  exit  pressure  ratio  is  plotted  as  a  solid  line. 
The  experimental  data  are  represented  by  the  individual 
points.  These  graphs  provide  a  sequential  look  at  what 
happens  to  the  flow  as  the  reservoir  pressure  is  increased. 
The  data  indicate  that  the  shock  in  the  reed  valve  occurs  at 
a  smaller  radial  distance  than  the  theory  predicts.  The 
daua  along  the  supersonic  portions  of  the  curve  also  match 
with  smaller  area  ratios.  This  appears  to  be  consistent 
with  the  idea  of  a  blockage  effect  in  the  flow  similar  to  a 
displacement  thickness  of  a  boundary  layer  (cl.  Chapter  11;. 
it  should  be  noted  that  the  flow  in  the  supersonic  section 
of  the  vaive  is  influenced  by  an  extreme  pressure  gradient. 
After  the  shock,  the  pressure  gradient  is  adverse  due  to  the 
suosonic  nature  of  the  problem  and  the  geometry  m  this 
region.  Also  of  interest  is  that  the  flow  in  tne  suosonic 
section  of  the  test  rig  actually  rose  above  ambient 
pressure.  The  theory  shows  that  the  pressure  ratio  should 
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Figure  20. 


Pressure  Ratio  vs.  Location  Plots,  Theory  and 
Data,  for  Reservoir  Pressure  of  15. QS  psia 
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Figure  22. 
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Figure  24. 


Pressure  Ratio  vs.  Location  Plots,  Theory  and 
Data,  for  Reservoir  Pressure  of  36. 39  psia 
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Figure  26. 


Pressure  Ratio  vs.  Location  Plots,  Theory  and 
Data,  for  Reservoir  Pressure  of  72.  90  psia 
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Figure  27. 


Pressure  Ratio  vs.  Location  Plots,  Theory  and 
Data,  for  Reservoir  Pressure  of  96.  433  psia 
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Figure  28. 


Pressure  Ratio  vs.  Location  Plots,  Theory  and 
Data,  for  Reservoir  Pressure  of  104. 9b  psia 


approacn  the  pressure  ratio  at  the  exit  from  beiow,  with  no 


overshoot.  The  data  appears  nearly  flat,  but  close 
inspection  reveals  that  the  measured  pressure  ratio 
approaches  the  exit  pressure  ratio  from  above. 

At  reservoir  pressures  above  170-inchea  of  mercury 
(pressure  ratio  of  approximately  0.15)  the  valve  plate  began 
vibrating  and  the  centering  pms  had  to  be  held  in  place. 
At  gauge  pressures  above  190- inches  of  mercury,  the  r lg 
emmitted  a  howling  noise  as  the  vibrations  increased. 
Pressures  above  210-inches  of  mercury  were  not  attempted  due 
to  the  intense  noise  and  instability. 

Figures  30A  and  30B  are  included  for  a  comparison  of 
the  data  collected  with  the  reed-valve  model  to  data 
collected  in  a  Laval  nozzle  (10:129).  A  striking  similarity 
is  apparent.  The  area  ratio  attained  by  the  reed  valve  is 
mucn  greater,  but  the  jump  in  pressure  due  to  a  shock  and 
t ne  overshoot  and  return  to  ambient  pressure  near  the  exit 
are  displayed  in  both.  Further,  it  should  oe  noted  tnat 
Figure  J0B  is  given  as  an  example  oi  how  well  the  one- 
dimensional  theory  holds  true  in  experiment  «10:129>. 

i'he  neignt  of  tne  vaive  plate  was  measured  over  a  large 
range  of  pressures.  The  heignt  was  measured  with  an 
automotive  feeler  gauge  and  found  to  be  constant  to  within 
the  tolerance  of  this  measurement  tecnmque.  The  heignt  was 
found  to  be  approximately  0.010-inches  for  the  range  of  tne 
test  runs.  The  bourdon  pressure  gauge  was  hooked  to  tne 


inlet  line  ol  the  tlow  in  order  to  make  a  static  pressure 


reading .  By  forming  the  ratio  of  static  inlet  pressure  to 
total  pressure,  the  Mach  number  of  the  flow  in  the  inlet  was 
found  to  be  approximately  0.  205-  Although  it  is  most  likely 
that  the  flow  in  the  inlet  is  turbulent,  a  quick  look  at  the 
problem  assuming  a  laminar  Hagen-Poiseuille  profile  (11:12) 
is  a  good  order  of  magnitude  check.  The  volumetric  flow  is 
given  by 

Q=(n/ 2J  f(,  but*  u,  (17) 

Where  u*  is  the  maximum  velocity  of  the  flow.  Recalling, 
that  continuity  may  be  written  as 

m  a  P  Q  =  p AV  <,  la  > 

The  throat  area  may  be  chosen  as  the  second  location  to  do 
the  comparison  where  A  =  2nRiNUETh.  Eq  (la)  is  then  solved 
for  the  height  h. 

h  —  (  p,  Ui  K,  phit  )  /  (4  p,  u,  )  1  a ) 

.-or  the  Mac.n  number  of  the  inlet,  the  calculated  heignt  is 
n  =  0 . 0 iWb - inches . 

-onxiguration  2 

The  1.5-inch  by  1.5-inch  square  piate  was  tested  ior 
gauge  pressures  ranging  from  a. 2-inches  of  mercury  to  140- 


mches  Ql 


mercury. 


The  plate  became  unstable  at  the  upper 


pressure.  The  plate  was  centered  above  the  imet  at  various 
orientations,  This  was  made  possible  by  drilling  additional 
alignment  holes  in  the  overhead  plate  through  which  the 
centering  pins  could  be  inserted  (See  Figure  12).  Pressure 
measurements  were  taken  with  the  square  aligned  with  the 
axes  of  the  radial  surface  ports  and  at  -15-degrees  and  at 
*J0-degrees.  The  pressure  profiles  were  expected  to  be 
different  for  different  radial  directions  beneath  the  non- 
circular  plates  due  to  the  variations  in  the  lengths  ox  the 
radial  rays  from  inlet  to  the  plate  edge. 

It  came  as  a  bit  of  a  surprise  that  the  data  appeared 
to  show  little  variation  with  respect  to  radial  direction. 
The  measurements  taken  with  the  plate  turned  either  -15- 
degrees  or  *30-degrees  also  showed  little  variation.  The 
data  has  the  same  appearance  as  the  data  gathered  lor  the 
circular  disk.  KEDDATA  was  put  to  use  to  determine  if 
indeed  there  were  variations  in  the  data.  Figure  31 
contains  three  pressure  profiles.  The  first  is  obtained  by 
averaging  data  from  the  0-degree  and  30  -  degree  radiaxs.  i’he 
second  profile  is  pressure  ratio  data  from  the  lour  45- 
degree  radials.  The  final  curve  is  plotted  as  n  tne  data 
were  Xrom  a  circular  disk;  all  eight  radial  profiles  are 
averaged  to  give  this  final  curve.  The  three  curves  are 
nearly  identical.  Obviously,  the  distance  travelled  in  the 
cnannel  tor  this  configuration  did  not  effect  the  behavior 
of  the  flow.  If  one  postulates  that  there  is  a  significant 


Pressure  Ratio  vs.  Location  Data  tor 
Conliguration  2:  Three  Data  Plots;  Two  tor 

Symmetric  Arms,  One  for  Circular  Evaluation. 
Theoretical  Plot  for  Circular  Plate  Included 
for  Comparison 


reparation  m  the  ilow  downstream  ot  the  shock,  then  the 

pressure  in  the  region  from  the  shock  to  the  exit  could  tie 

close  to  ambient  along  all  ray  directions.  The  pressure 

ratio  curves  are  quite  flat  in  this  region  so  this  seems 

like  a  reasonable  possibility  of  what  is  happening.  If  this 
is  the  case,  then  a  longer  or  shorter  channel  in  this 
portion  of  the  flow  will  probably  have  little  effect  on  the 

overall  characteristic  of  the  flow  in  general.  The  geometry 

of  the  inlet  is  circular,  the  supersonic  flow  issues 

radially  outward  until  it  passes  through  what  appears  to  be 
a  circular  jump  discontinuity.  It  seems  likely  that  the 
geometry  of  the  inlet  rather  than  the  geometry  of  the  valve 
plate  disk  is  the  determining  factor  of  the  ilow. 

It  should  be  noted  that  the  data  on  the  supersonic 
portion  oi  the  curve  for  the  square  plate  ao  not  tall  on  the 
same  curve  as  the  points  for  the  circular  plate.  The 

ei  feet  i  ve  length  of  the  flow  cnanr.ei  in  eacn  case  is 
diiierent.  The  points  are  close,  but  a  distinct  difference 
is  evident.  The  theorectical - isentropic  subsonic  to 

supersonic  portion  tor  tne  one-inch  diameter  circular  plate 
is  included  in  t-'iqure  Jl  lor  comparison  oniy. 

configuration  o 

This  rectangular  valve  plate  measured  1. 75-inches  by  1- 
inch.  Test  runs  were  performed  for  gauge  pressures  from 
d- inches  ot  mercury  to  ia«-inches  of  mercury.  The  plate 
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PRESSURE  RATIO  VS  LOCATION 


RADIAL  DISTANCE  (INCHES) 


Figure  33.  Pressure  Ratio  vs.  Location  Plots  of  Data  and 

Three  Analytic  Techniques:  (1)  SOLVE, 

<2)  EFFAREA  Method  1,  (3)  EFFAREA  Method  2. 

Reservoir  Pressure  is  36. 39  psia 


larger  radiai  distance. 


Method  1  does  match  the  data  points 


on  the  supersonic  curve  fairly  well.  But,  this  should  be  no 
surprise,  because  the  area  ratios  were  derived  from  the 
supersonic  portions  of  the  data  plots. 

Method  2  predicts  the  shock  location  well  using  a  value 
of  to4. 0  percent  for  the  effective  area  of  the  throat.  This 
value  for  the  effective  area  matches  the  shock  fairly  well 
for  the  r ange  of  pressures  tested .  However,  the  values 
along  the  supersonic  portion  of  the  curve  are  incorrect .  In 
fact,  these  values  are  worse  than  those  predicted  by  SOLVE. 

Neither  method  of  EFFAREA  is  adequate  to  predict  the 
flow.  More  work  to  attempt  to  match  the  aata  is  probably  m 
order,  perhaps  in  a  future  study. 


Flow  Visualisation 

The  oil  now-visualization  procedure  gave  a  distinct 
depiction  of  where  the  shock  was  located.  A  line  was 
produced  in  the  oil  by  the  shocK  in  most  situations. 
Painting  the  base  plate  of  the  test  ng  black  roughened  the 
surface  and  slightly  reduced  the  sharpness  of  the  shock 
line. 

For  coni iguration  1,  three  test  pressures  were  chosen, 
the  test  pressures  were  44,  73,  and  33  psia.  These 
pressures  were  chosen  from  a  study  on  the  reduced  pressure 
ratio  data  already  gathered.  The  shock  locations  should  be 

i, 


between  ports  1  and 


and  3,  and  3  and  4  respectively. 


Figures  34A  and  J4B  show  the  flow  visualization  ior  tne 
first  two  runs.  Each  shock  was  well  defined  and  at  the 
expected  location.  The  shock  in  each  case  was  circular. 
Figure  35  shows  the  oil  flow  for  the  99-psia  run.  This 
pressure  corresponded  to  about  170-inches  of  mercury  and,  as 


reported  earlier,  the  plate  began  vibrating.  Close 
inspection  of  the  oil  flow  does  reveal  that  the  shock  was 
located  at  or  just  beyond  port  3.  However,  the  on 
visualization  is  smeared  as  would  be  expected  if  the  shocn 
were  oscillating. 

One  test  run  was  performed  with  configuration  2.  The 
oil  flow  is  shown  m  Figure  36.  The  test  pressure  was  39- 
psia.  Although  the  plate  is  square,  the  shock  was  circular. 
This  supports  the  idea  of  separation  discussed  earlier 


(cx.  Chapter  VI  Coni iguration  2>. 

Two  test  runs  were  performed  with  the  rectangular  plate 
configuration  3).  The  first  run  was  for  a  reservoir 

pressure  of  39-psia  and  the  flow  visualization  is  presented 
m  figure  37A.  Again,  the  shock  was  circular.  The  second 
run  was  pertormed  to  place  tne  shock  near  the  iiose  edge  ox 

the  prate.  A  pressure  was  chosen  just  below  the  threshold 

wnere  tne  pra;t  began  to  vibrate.  The  reservoir  pressure 
was  o4-psia.  The  snocK  was  again  circular  as  seen  in 


Figure  j  7 b . 
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The  on  drop  flow  visualisation  worked  well  but  was  not 
suited  tor  looking  at  tne  shock  alter  tne  vibrations  set  in- 
Possibly  a  Sctilieren  technique  with  a  plexiglass  model  would 
be  a  test  bed  to  look  at  the  higher  pressure  flows. 


Discussion 

The  one-dimensional  isentropic,  inviscio  model  with 
normal  shocks  imbedded  in  the  flow  predicted  the  basic 
nature  ot  the  flow;  in  tact,  as  well  as  they  predict  lor  tne 
uavai  nozzle  in  Reference  lid. 

it  is  interesting  that  this  reed-vaive  mcdel  so  closely 
matched  the  performance  of  a  Laval  nozzle  in  the  supersonic 
section,  especially  considering  the  equivalent  divergence 
angle  oi  the  reed-valve  model  (cf.  Chapter  IX).  Hill  ana 
rsrerson  i  4  : -I'd  / -  40  /  >  state  that  the  best  performance  is 
obtained  tor  conical  nozzles  with  divergence  naif  ar.gies 
between  12  and  lb-degrees.  According  to  Reference  4,  this 
preferred  range  of  half  angles  is  due  to  tne  stated  "fact" 
that  larger  angles  nave  the  possibility  ot  now  separation 
from  the  wans. 


vii.  conclusions  and  Recommendations 

conclusions 

The  purpose  of  this  project  was  to  analytically  predict 
and  experimentally  measure  the  flow-induced  pressure  field 
in  a  reed-valve  arrangement.  This  pressure  field  was 
initially  expected  to  be  that  lor  fully  subsonic  flow.  The 
flow  in  both  the  reed-valve  analytic  model  and  experimental 
reed  valve  was  supersonic  in  nature  with  imDedoeo  shocks  for 
all  but  the  lowest  total-pressure  runs.  Even  the  spool  and 
paper  trick  can  readily  be  made  to  have  supersonic  flow,  a 
normal  shock  and  finally  a  subsonic  flow  to  the  exit. 

The  one-dimensional  isentropic-f low  model  with  an 
imoedded  normal  shock  is  a  good  model  for  this  flow.  There 
is  a  definite  phenomena  in  the  flow  channel  which  looks  like 
a  biocxage  due  to  tne  displacement  thicxness  of  a  Dounaary 
layer,  separation,  or  both.  The  phenomena  did  not  yield  to 
simple  modeling  techniques.  The  addition  of  a  model  to 
account  for  boundary  layer  growth  might  further  enhance  tne 
model . 

The  effective  area  methods  of  analyzing  tne  problem 
need  more  refinement.  The  task  of  closing  the  reed-valve 
problem  to  predict  the  valve  height  and  tne  mass  flow  rate 
can  be  approached  from  simplified  models  to  get  a  fair  order 
of  magnitude  approximation. 

The  flow  visualization  technique  worxed  well  lor 
getting  a  definite  shocx  location  and  shape.  Krom  the 


smearing  oi  the  shock 


line  lar 


test  runs 


when  t ne 
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piate  was  vioratmq,  it  seems  that  tne  vibration 
an  osci^iation  oi  the  shock.  A  Sch neren  tecnm 
effective  for  determining  the  shock  character isti 
high  pressure  flows. 
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The  analysis  o±  the  problem  13  certainly 
with  a  small  computer  system  when  using  the 
technique. 


Recommendations 

Based  on  the  experiences  and  the  results  01 
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!  tne  iol.oving  recommendations  are 

i.  investigate  methods  to 

tne  problem  to  determine  the  neiqh 
tne  case  plate.  A  momentum  approa 
_.  investigate  bounaary-i 

j033iDi.it y  ox  separation  m  the  x 
predicts  the  effective  area  ox 
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considerations  in  mind. 


Predict  mass-tiow  rates  from  the  above  extensions 


j . 

to  tr.e  modelling.  Verily  the  procedure  by  experimentally 
measuring  the  mass  tlow  rate. 

4.  investigate  various  inlet  geometries  to  find  the 
eitect  on  the  supersonic  flow  and  the  shape  of  the  shocK 
wave. 


b.  Use  a  Schlieren  system  in  conjuction  with  a  clear 


plexiglass  reed-valve  model  to  more  precisely  define  tne 
character  of  the  shocK  at  higher  pressures. 
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Appendix  a 


eueVE 

A  complete  listing  oi  the  program  eOeVE  lollows.  sulVl 
analytically  solves  the  reed-valve  problem  lor  a  circular 
valve  plate  and  circular  inlet.  The  assumptions  wmch  SOLVE 
uses  are  noted  m  Chapter  II. 


1000  ’  PROGRAM  -  SULVE 

1010  ' 

i0.-0  '  T  HIS  PRGljKAM  ANAL  /  TIlALl  <’  bULVha  THe  ri_  — w  r  ;lLj 

I'ilbU  '  UNDER  A  CIRCULAR  L)lbK  USING  1  -  D  ibLNTKUr'it 
1040  'AND  NuRMAL  SHUCK  FLOW  RELATIONS. 

.000  ’ 

1060  UiM  R< 19 > ,  P < 19) ,  M( IS) ,  V( 19) , C < IS ) ,  RHUi  J S/ ,  Tt  IS > .  PR i  Is > 
1070  ' 

10U0  '  DEFINITIONS 

10S0  ' 

1100  GAMM A  =  1.4:'  RATIO  OF  SPECIFIC  HEATS 

1110  RO  =  1714.7:'  GAS  CONST  iFT''2/DEG  RANKlNE 

1110  PA= 14. 695b : '  AMBIENT  PRESSURE  iPSI) 

iij0  TA  =  iiiS:  '  AMBIENT  TEMP  iRANKlNE, 

-  -  *40  I  i  =  i  A  :  '  TOTAL  TEMPERATURE 

i.D0  R 1 -  .  Led:  '  INLET  RADIUS  l  INCHES ) 

i  ibid 
i.i/0 

11S0  '  INPUT  CONDITIONS 

i  -  3ld 

1-00  PRiNi  "DEFAULT  UPTiU.NS  ARE  FOR  A  STD  DAY 
1-10  iNPU  T"DU  YOU  WAN  T  TO  USE  DEFAULT  UP  T iUNs  T  uls 

*  i*  J  1  s  =  "  N  "  uul  iO  i-Aw 

1 . _ 3  0  ip'  iJii="n"  uuTu  l-NW  sLsE  uu  Tu  l_/0 

-  _*40 


J. 

ijusub  -li0 

:  '  GO  TO 

ALTERNATE  INPUT  bUbKLiUi 

i'iNt 

* 

A.  -  /  W 

iNr'J  .  *  iNPU 

T  ;  HE  KAiJiJb  uK 

THE  uibft  ii.NcrtEc.1  :  ”, 

n  E 

i  -au 

.  NPU T* 1 NPU 

T  THE  i N Lc.  T  PkEb 

SURE  i N  FblA  : 

PRES 

- 

r-  .  r  r  tt  E  S 

^  -  Id'  u 

i  -3  i  wj 

uEh : NE 

fune r iuNb 

1 

Dc.r  r  NbuBul  )  =  (  i  ♦  (  i  bAMMA-  1  ) 

/  2  ;  *  M  ■'  ..  > 

' 

i. 

DEFINE 

EXPONENTS 

1 

e X 1 -DAMMA/ 

(  UAMMA-  1  i 

.-■bid  s  X  _  =  1  /  i  uAl'iflA-i  ; 


-  s  /  0  s  A  -j  =  .  O  -  i  UAtir.A  /  (UAllMA-i.  )  > 


1260  EX4  = i OAMM A  -  1 ) /GAMMA 
1260  * 

1a00  *  utFiNii  TOTAL  PARAMETERS 

1410  c  i’  =  l  GAMMA  *  RG*  l’T  )  ■' .  5  :  '  SUNiL  V  EL  FOR  TuTAL  lUNU 

*420  RHUT  =  PT/ ( RG*TT )  :  '  TOTAL  DENSITY 

1400  ' 

14443  '  INCREMENTAL  LOCATIONS  FOR  CALCULATIONS 
1450  R= < ( RE-RI ) *16) -1 

1460  I NCR= 1 / 16 : '  EVERY  1/16  INCH 

1470  FOR  1=1  TO  R 

1460  R(I)=Ri*(I-l>* I NCR 

1460  NEXT  I 

1 300  ' 

1510  '  DETERMINE  PT/PE  FOR  SONIC  THROAT 
1520  ' 

1520  '  FIND  MES  -  MACH  (3  EXIT  FOR  j'!’A  R  lOND 

1.-340  MEs  =  .  i  :  '  JUST  A  FIRST  HACK  To  STAR  T  ITERATING 

1220  MES  1  =  <  R I  /  RE  »  (FNBOB(  MES)  /FN  BOB  l  1  )  )  •'  (  -  EX  2  ; 

1260  IF  AUS iMES-MEsl J <. 000001  THEN  GOTO  1560 
*2/0  MEs=MEsi 
1260  GOTO  1550 
1560  MES  =  MEs 1 

1600  '  NOW  USE  MES  TO  FIND  PT/PE  FOR  STAR  LUND 

.1610  '  CALL  IT  PRS  —  FOR  PRESSURE  RATIO  STAk 

1620  PRS= i FNBGB<  MES)  ) A  EX  1 
lb  20  ' 

1640  '  CuMPARE  PRESSURE  RATIO,  PR,  WITH  PRS 


1620 

r1  K  - 

/  PA  :  ' 

PR  Is  PRESSURE 

R  A  T I  O 

1660 

r  K 

>  PKS 

THEN  OUTU  1660 

ELSE  DUTO  1/20 

1  6  /0 

' 

i  660 

UUij'J  o 

2*120  :  ' 

GUTO  1’HE 

sUPERSONic  sUbruuTiNe 

1  660 

jU  I  U 

6620 : ' 

GUTO  T 

HE 

OUTPUT’  SECTION 

1  ‘00  ' 

*  .'  10  ' 

*/*0  '  jUtlbUN  1C  SOLUTION 

.  /  _  0  r  .4  -  N  . 

*  /  “  0  R  i  N  I  ’’  TH  1 6  i  S  A  sUBsuNic  PROBLEM 
.  /  jiii  P  R  1  N  1 

*  60 

1  /  /0  .-.E=  •  '2  /  ••  GAMMA  -  1  ;  )  *  4  v  PR  A  EX4  )  -  I  )  .'  A  .  a  :  ‘  scHREiER  146 

L  60 

*  /  60  Wii’H  ME  IN  BAND  FiND  :1  ;  i  i  Ai  All  kAuiAl  .ulauuN 

1 600  'JK  *  ~  i  .  U  ri 

1  6  1  0  LjUSUB  1  6  1  0 

i6_0  Nee  i  l 

ia_>0 

1640  WITH  Mi's  ut T  ace  uihEn  PRUPeRTIes 

.620  rUK  1  =  1  i'J  ei 
1660  UUsUB  2020 

1670  NEXT  1 
1660  UO  TO  2620 :  ' 


GOTO  THE  OUTPUT  SEcTiUN 


x  SOW 

IslG 

i 

i-scG 
iS4G 
isSG 
1  goO 
ly/u 
lyaw 
lsso 
uoiDC 
S01G 
cUx.g 
x  WJ  _>u) 
CG4G 
Oobtd 

-.LJSU 
-.Id  /id 
OdtUd 
-  GSld 
Gild  Id 
Cl  1(J 
x-icld 

■iiJW 

C140 


S UBHUUTIUE  --  INCREMENTAL  MACH  NUMBER 
INPUT  --  EXIT  MACri  NUMrisR 
fli  =  .  3  :  1  ANuTHEk  HACK  TO  START  THE  x'i  ERA  i  luN 
)  *  (  P’NBGB  (  ME  )  / FNBOB  (fll)  >  AEX 0 
IE  AUSi  Ml  -  Mil  i  < .  ididididid  1  THEN  GuTG  lyaiJ 
MI  =  M 1 1 
liU  l  U  lS4ld 
M( I ) =  M 1 1 
RETURN 


SUBROUTINE  --  INCREMENTAL  PROPERTIES 


INPUT  -- 
BUB=FNHUB  (  M ( I >  > 
p  i  I > =PT/ l BOB  ' EX 1  -  :  ’ 
rK i : ) =Pl I ) / PREs : 

KHO i. i > =RHUT/ ( BOB  A  EX  I 
AnOt i i =HHUt  I  .■  *  lidOG  :  ’ 
f i 1 ; - TT / BOB : ' 

Ci  x  > =  O T /  i BOB  A . i )  :  ' 

V l i i =Mi I j  *Oi  I  i  :  ' 
RETURN 


MACH  NUMBER 


PRESSURE  iPSiA) 

PStaaliKt  «A  iTu 

DENS  IT?  i SLUGS/ cUBic  Eu 

RriUC  X  f  E  - 'd  s 

TEMP  i DEGREES  RAfiAiNt  j 
SONIC  VELOCITY  C  FT/SEC) 
FLUID  VELOCITY  i FT/ SEC ) 


Euc  1  > 


-jlbtf 

' 

SUBROUTINE  --  ALTERNATE  INPUT 

INPUT " 

INPUT 

GAMMA  : 

^  i  /  VJ 

x  N  P  U  i  " 

i  NPUT 

THE  TuTAl  TEMP  IN 

uEGKt.Es 

a  : 

ciCW 

.  N  P  U  I  " 

xNpu  r 

the  exit  pressure 

IN  P  s 1 A 

: 

*  MrU T " 

INPUT 

THE  R  A L> 1 U  S  G  F  THE 

INLET  C 

I  N  C  H  c. 

RE  TURN 

'  ,  GAMMA 


sUBROUTxNt  --  sUp ERScN  x  c  scx-c  i  xUN 


r  R  x  N  T 

-x  X_  '~J  kj 

p  r  .  n  ;  " 

This  x  s  A  sUPSRScNic  pi- 

lubes/. 

— 

P  R  i  NT 

/kJ 

3  U  P  s  R  B 

- ” ILg” 

.  -  *.  o  «J 

G  1  R  =  X  : 

-Nx  llAx  MArcH  x'j 

lRSCixuN  =  xN 

P  Acs  =  . 

tu  :  ' 

g  i  c  E  Or  MARCH  s  ■ 

,'EP 

;  -  :x  *  N  .  •' 

i'IAaL  rl 

THE  oriLJ.rv  xN/CUT  UN  l  ie 

p  E  =  p  A  "  :  p  R  x 

1-  -  A.  W 

S  A  x  ,  . 

wJ kjvj  x  : 

PRs.ssUrs  u  Ip  FErsN  T  x  Ac  ?  sa  i  i 

A.-J  —W 

ft  A  r  : 

xNI  i’xAc  uUsss  ? 

fl  A  A  1N1  Aun 

r  r  -  ■  R  X 

/  M  AIM 

i  Il-AiIHA*  i  1  /  x  1  /  r  NBUBeIX 

t  >  '  *  C.  A  J 

_  4id 

_ '  ID  Vj 

x  X  /  X) 

_  :SCI 
..cGId 
-  i  ixi  w 


fir  •  i  i  i  uAfIMA-i  /  ‘flX'x'x  W  l  *bAMMA*flA  T  )  -  i  bAilllA-H  i  i  '  . 
uluBI  =  (  i  i  ■  UAilM  A  *  i  >  /  G  )  *  MX  '  0  ;  /  FNBOH  ■.  MX  )  )  '  EX1 
jluBx=i  x  GAMMA  *  it/  X  x  c*GAM.1A»Ma  I  -  \  uAMMA-  1  l  >  >  '  s  A  C 
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—  440  Me.=  MEi 

2450  uu  :  U  -420 

24o0  ME=ME1 

24/0  ' 

4ts0  PE=PTY/ i FNBOBc ME ) ) AEX1 
4  iw  ' 

2500  '  COMPARE  EXIT  PRESSURE  WITH  AMBIENT 

2510  '  MCJVE  SHOCK  UNTIL  THEY  MATCH 

2520  IF  ABS i PE-PA ) <EXIT  THEN  GOTO  2750 

25J0  iF  PE<PA  THEN  GOTO  2650 

2540  ' 

2550  '  MARCH  OUT  WITH  THE  SHUlK 

2560  MARCH  =  0:PRINT"MARCHING  OUT  --  PE/PA  =  ",  PE/PA 
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-  RElAl-ULA 

-O00 

palE  =  paCE*  »  .  5; 

2bR  =  0 
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2050  '  MARCH  IN  WITH  THE  SHOCK 

26tr.0  MARoH=I  : PRINT "  MARCHING  IN - PE/PA  =  ".PE/PA 

20/0  IF  MaRlH-UJ.R=0  THEN  2660  ELSE  2700 
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4000  I p  aUFEki  =  "  YES "  THEN  4010  ELbE  40b0 

40iU  PRINT  «V ,  tJKX 

4030  PRINT  RY,PRY 

4030  WRITE  #1, RY, PRX 

4040  WRITE  #1, RY, PRY 

4030  ' 

4060  FOR  I  =  ISH<-1  TO  R 
4070  PRINT  R  (  I  )  , PR ( I ) 

40b0  WRITE  #1,  R l I > , PR ( I  ) 

40b0  NEXT  I 
4100  CL.ObE 
41 J0  PRINT 

4100  iNPU'1  “  bHAU.  I  SAVE  i  HE  pRESbURE  uAfA  "  ,  UJ'l 
4 1 00  IF  yjS= "y*  THEN  4ioU  ELbE  4 140 
•JitW  .  F  Wjb  =  ”Y"  THEN  4io0  ELSn.  3000 
‘UjiJ 

-it.0  '  UPEN  l»A  i'A  FILE  FUR  PRESSURE  Vb  K  DATA 

4170  ’ 

4160  INPUT'INPUT  THE  NAME  OF  THE  DATA  FILE  :  M,Li0S> 
4 1  *40  PRINT 

4000  OPEN  "U", #1, u2S 
4-iU  '  WRITE:  R( I > ,  r ( I > 

4000  WRITE  #1,0, PRES 
4 330  PRINT 

4-40  P  K 1 N 1  "  PhtbbURE  DATA 

4-30  PR. NT 

4000  pkiNT"  k 1 1 :  Ptl> 

4-/0  print  0, pres 
4-60  fur  1=1  Tu  10H 

4-bU  PRINT  K  1  i  )  ,  P  i,  :  ) 
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Appendix  b 


instrumentation  List. 

item  serial  Number 

electronic  Control  Tecfinology  S-100  Computer  - 

Anaiog-to-Digital  Card,  Dual  Systems  Control 

■loaei  Aiflli  REV  D .  ^.-SjlS 

jiaitai-to-Analog  Cara,  Dual  Systems  Control 

rtociex  AuM-ll . j.uja 

i.csn  Jnve  Assembly,  Tarbell .  ilb  /  - i<U4a.L 

iermmai,  He  a  true  it  fioaei  H-19 . .  1 - 4  j V  i2 

team  valve,  SoANuu  Model  4BS9  3ti)12 . 

Transducer,  SCANCQ  Model  PDCR23D- 25PSID  .  .  .  - 

uuai  uc  Power  Supply,  Hewlett  Packard  s2tabB  .  il4wAtd9w40 

Solenoid  Controller,  SCANCQ  CTLR2/32-SS  .  .  .  14a2 

odd-tven  uecooer,  slANlu  uEU^. .  /'«*odcu 

iWC-mcn  Hg  Manometer,  King  engineering  .  .  .  w*.Oz 

.WB-:ncn  rtg  equivalent  bourdon- type 

Pressure  uauge,  Wallace  and  Tiernan . .-A  -  40 - nn x wuue 

Jiaita*  Multimeter,  mew^ett  PacKara 
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TEST 

A  complete  listing  oi  program  T£bT  tonows.  i  ms 
program  was  used  lor  data  acquisition.  A  description  ol  the 
application  ol  TEST  is  given  in  Chapter  XV. 
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i  7  id 
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-1  <d 
_ C 

—  fciU 
_  /  iJ 
_ t5  ud 

_  1  'J 
-5— 
_5t*j 

-0*4 
/  4' 
JuU 

-J 

‘ivd?«u 
4  iW 
H-'O 
■4  ju 
440 
400 
4cjC 


'  PROGRAM  --  TEaT 

THIS  PROGRAM  GATHERS  DATA  USING  A  SCANIVALv'E. 
i  I  3  i'h.r’i  J’HE  VaevE  i'HRuUuH  PORTS  ANu  3Ai'ir»c.3 

EACH  PURI'  "h"  i'iilES.  1'HE  AVERAGE  OF  "HE  "  K " 
SAMPLINGS  ;3  PRINTED  FOR  EACH  ruR  i  . 

'  bt-P'GRE  RUNNING  THIS  PROGRAM  FOR  l>A1  A.  uR.E-n 
THE  POWER  SUPPLY  OUTPUT  AND  THE  CORRESPONDING 
OFFSET  v  IOV  AND  00WS  ). 

E  =  LO:'  #  OF  SAMPLES  TO  BE  AVERAGED 

INPUT” INPUT  CONFIGURATION  #  :  " , CONFIG 

xNPUT" INPUT  RESERVOIR  PRESSURE  :  ”,PU 
ro  =  RO~  .  i 
R  R  i  N  Y 

R  siNi 

rRiN  i'  "uUNr  IGURA  I’iUN" ,  CONFIG,  "RES  PRtaa  “in  :ig 

R  R  i  n  i 

fxINT"  rOR  i's  i  I'HruUljH  i 

PRINT"  RADIAL  VOLTAGE  IN  JIIll.  'VoL  Is 

R  K  i  N  'l 

For  .kau  =1  10  a 

r  un  .  Poh  i  -  i  I'O  3 


uu  r  i  /  /  , 

i  *.  /  :  ' 

VOL  r AGE 

u  N  —  v  a  _  v  b 

uUT  X / ^ » 

1x7 

l  =  d:  ' 

JtLA r  r > 

j  A  co  *  x  .  i_  n  ^  «\  o 

r  Jh  u  j.  =  ^ 

.  U  L 

NLA  .  x  i 

uU  ;  ^  f  , 

J  :  ' 

'/Uu  i'Abt 

ur  r 

uU  .  x  b  , 

G 

r’jn  A  -  ~  i 

|U  L. 

imex  r  a  i 

oL.':. 

goos 

A  A  U  U  =  G 

FUR  1  =1  I'O  K 


<iNP(14u)  ANL>  ID)  »13b 


-te>0  uli  i  i4t),  i- hi 
**U0  UU  r  l-il  ,  0 
3W0  A  =  INF  i  14*:  )  » 

:Di0  A  =  A  -  20413 
3*.0  rt=  ;  1  0/  40SS  )  *  A/  100 
D5U0  /»ADD  =  AaDD*A 
540  NEXT  I 
550  AAVE=AADD/K 

it>0  A  AD  J  (  ikAD,  IPORT  )  =  (  AAVE-UFFSET  )  »  1000 
570  NEXT  IPGRT 
■300  PRINT  USING 
;ei  PRINT  USING 
5a4  PRINT  USING 
D.UU  NEXT  ikAD 

O  00  ' 

c 10  - ur  i-= i  ru 
oj-0  JU  T  1  7  /  ,  127 
o.:0  jU  r  i/b,  127 

o4u  .-UK  1  *:  =  1  lu 

eD0  NEXT  J2 
eo0  ULI  17  7,  0 
o  70  uU I  i 7b, 0 
bb0  FUR  K2=l  Tu 
S90  NEXT  K2 
700  NEXT  12 
710  * 

720  '  SAVE  DATA  TO  DATA  FILE 

■  1 0  F  k  i  N  i 

/‘to  :  N  F  U  !  "bo  I'UU  WANT  fU  oAVE  TRIf  uh  T  A  •'  ”,  ui 
20  -1  0i=  "V"  UU  i'U  O00 

7»u  iNPUi” HOW  '  BuU  T  ANOTHER  RUN  ",  Uli 
7/0  _ r  0*.s>=  "  I  "  ulilU  100 

cU  fRj.Ni  "  i  i-  '  Sa  bEtN  A  FLEAF.URE  DUiNu  bUF.iNe.oo  Wxi'H  llj 
:  '0  UO'i'Li  ;  lt) 

300  fain  !  "  INPUT  FILE  NAnt :  d  :  Rflfluu#* .  l>A  . 

.3*0  *  ,h  i  u  i  "  .-  *  j.  a.  n  A  il  E  :  " ,  i- 1> 

■13  *.  0  -J  F  e.  N  "U",  H,  ra 

a u0  wRi.Eff  1,  ouhr'.u 

0*3  0  W  R  *  i  i»  *  ,  F0 
_'7j0  r  UK  j  u  -  i  .U  b 
dab  r  UK  1j-i  i'U  U 
■3.0  HRiill  i,  AAl.  v  ja,  *Jl 
:o  U  .teA  i'  *  u 
o'-a  Ne.  A  .  J  i 

F  0  0  —j-UF.il  * 
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rEDDAT  A 

A  complete  listing  ot  program  kEDOATa  xoilows.  mis 
program  was  used  to  reduce  the  data  which  were  generated  by 
program  TEST.  The  details  ol  tne  program  are  described  in 
Chapter  V. 


LUOO  PRUuRAM  —  REDDATA 

iaiu  rrINT 

*0^0  r'KiH  f M  f'hL'UrtAfl  AtwEssita  DA  A  r'lLEb  btNt.tA  t'c.u  a  i  .  u  i  ' 
Idiiii  PRINT "THE  DA  i'A  ARE  PRESENTED  iN  AUSjLU'tE  ."KiasUKLa 

aLuPa  =  . 7Cbi SLOPE  Of  CAL  uuKVt  i  mV  Tu  in.  ng ; 

iUbU  PRINT 

.td/'O  INPUT  "ENTER  .DAT  FILE  NAME  :  » , GS 

J.UG0  PRINT 
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